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ABSTRACT 

We present a sample of 1899 galaxies with a close companion taken from the Sloan Digital 
Sky Survey Data Release 7. The galaxy pairs are selected to have velocity differences Av< 300 
km s" 1 , projected separations (r p ) < 80 hzl kpc, mass ratios between 0.1 and 10, and robust 
measurements of star formation rates and gas-phase metallicities. We match the galaxies in 
total stellar mass, redshift, and local density to a set of 10 control galaxies per pair galaxy. 
For each pair galaxy we can therefore calculate the statistical change in star formation rate 
(SFR) and metallicity associated with the interaction process. Relative to the control sample, 
we find that galaxies in pairs show typical SFR enhancements that are, on average, 60% higher 
than the control sample at r p < 30 }cl kpc. It is at these small separations that the strongest 
enhancements in SFR (by up to a factor ~10) are measured, although such starbursts are 
rare, even amongst the closest pairs. In addition, the pairs demonstrate more modest SFR 
enhancements of ~30% out to at least 80 fcl kpc (the widest separations in our sample). This 
is the first time that enhanced SFRs have been robustly detected out to such large projected 
separations. Galaxies in both major and minor mergers show significant SFR enhancements 
at all r p , although the strongest starbursts (with SFR enhancements of a factor of ~ 10) appear 
to be found only in the major mergers. We also find evidence that SFR enhancements are 
synchronised in an interacting pair, such that a higher SFR in one galaxy is accompanied by 
an increased SFR in its companion. For the first time, we are also able to trace the metallicity 
changes in galaxy pairs as a function of projected separation. The metallicity is generally 
diluted in galaxy pairs by ~ 0.02 dex, with an average metallicity decrement of -0.03 dex at 
the smallest separations, a trend that mirrors the SFR enhancements as a function of r p . The 
SFR and metallicity trends with projected separation are interpreted through a comparison 
with theoretical models. These simulations indicate that the peak in SFR enhancements at 
small separations is due to systems near the end of the merger process. The extended plateau 
in SFR enhancements out to at least 80 Kjl kpc is dominated by galaxies that have made 
a pericentric passage and are now experiencing triggered star formation on their trajectory 
towards apogalacticon, or on a subsequent close approach. 
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1 INTRODUCTION 

Galaxies that experience a close encounter with a companion are 
expected to undergo significant changes. Both observations and 
simulations have been used to probe the internal properties of 
galaxies as they progress through a merger, with simulations pro- 
viding a framework through which the observational results may be 
interpreted. Currently, theoretical models present a consistent gen- 
eral picture of the evolution of a galaxy in a merger. Strong tidal 
interactions may trigger bar instabilities in the central regions of 
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the galaxy in both the stellar and gaseous components. These bars 
are misaligned, and torques exerted on the gas by the stars result in 
the loss of angular momentum in the gas at larger radii. This gas 
then falls towards the nucleus of the galaxy, efficiently funnelled 
by the bar instabilities (Barnes & Hernquist 1996, Mihos & Hern 



|quist|1996[|Cox et al.[2006; Di Matteo et al. 2007 , Montuo ri et al. 



|2010||Rupke et al.|2010a||Torrey et al.|2012) . Indirect signatures 
of this picture are visible in the gas-phase metallicities and star for- 
mation rates of the interacting pairs. Gas inflowing from the outer 
regions is generally of a lower metallicity than the nuclear regions, 
so inflow from larger radii results both in a diluted nuclear metallic- 
ity < |Kewley et al.|2006||Eilison et al.|2008b||Michel-Dansac et al.| 
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|2008[ > and in a flattening of the standard metallicity gradient present 
in spiral galaxies iRu pke et al.|2 010b, Ke wley et al.||2010| |Perez| 
[et~aI1|20TT) . This new central concentration of gas provides the 
ideal catalyst for a significant starburst. Higher than average central 
star formation rates occur in nearly all simulations of galaxy merg- 
ers (e.g.,|Mihos & Hernquist|[l9^4l [7996) |Di Matteo et al.|20(J7] 
|Montuori et al.|2010[>, and are also a ubiquito us feature in observa- 
tional studies (e.g.,|Larson & Tinsley|1978||Donzelli & Pastoriza| 
79971 [Barton et al.|2000HLambas et al.|2003UAlonso et al.|2004| 



2006; Woods & Geller 2007 ; Ellison et al. 2008b 2010 Darg et al. 
20Tol|Xu et al |2010) . 



The strength of the induced star formation can vary dramat- 
ically from merger to merger. A significant contributing factor to 
this variation is predicted to be the mass ratio of the two galax- 
ies ( |Cox et al.|2006[ >. Galaxies in major mergers, i.e., with masses 
within a factor of 3, have previously been reported to show the 
strongest star formation rate (SFR) enhancements, on average, with 
both galaxies showing comparable enhancements (Woods & Geller 
|2007||Ellison et al.|2008b"l|Xu et al.|2010||Lambas et al.|2012) . Ad- 
ditionally, simulations suggest that a number of parameters within 
the merger exert a significant influence upon the extent to which 
gas is funnelled to the central regions of a galaxy. Most notably, 
orbital parameters can control the strength of the SFR response ( |Di| 
|Matteo et al.|2007|[D 7 Onghia et al.|2010| >. The difference between 
a prograde-prograde and prograde-retrograde merger can alter the 
triggered star formation by a factor of 2, and the separations of the 
nuclei at first passage can inhibit SFR at coalescence if the tidal 
forces are so strong that gas is preferentially ejected into tidal fea- 
tures, rather than collecting in the central regions. Gas fractions 
also seem to play a role, with lower gas fraction galaxies showing 
weaker SFR triggering (Di Mat teo et al.|2007| >. 

Observationally, studies of the SFR in galaxy pairs with pro- 
jected separation (r p ) are largely in agreement with each other. 
Galaxies in pairs show enhanced SFRs out to separations of ~ 30 
kpc, with the strongest enhancement at the smallest separa- 



tions, and a smooth decline to a fiducial value (Lambas et al. 2003 
Alonso et al.|2006l|Nikolic et al.|20041|Ellison et al. 2008b; Li etal.| 



2008 ). Some studies have found evidence for SFR enhancements at 



wider separations (e.g., 


Barton et al. 2000||Lin et al.|2007||Robaina 


|et al.|2009; Wong etal. 


2011b, although generally this only weakly 



extends the trend to ~50 kpc. Galaxies are generally expected 
to be most strongly enhanced just after first passage, and again at 
coalescence, but the timescales for catching a galaxy during a sig- 
nificant starburst are also tied to the strength of the burst (Torrey| 
|et al.| [2012 l. The highest SFR enhancements are generally of the 
shortest duration, whereas lower-level enhancements can persist for 
longer periods of time l {Di Matteo et al.|2007{|Montuori et al.|20i0) . 
The trend of SFR with projected separation is therefore likely to 
be a combination of merger-induced changes and the timescale on 
which the observational snapshot is taken. 

In practice, our ability to determine the merger phase of a 
galaxy pair is severely limited, as our only observable indicators 
are that of the projected separations (r p ) of the two galaxies, or the 
morphological disturbance of the galaxy. However, even this latter 
metric is fraught with problems. Depending on the orbit of the en- 
counter, galaxies do not always show strong tidal features after an 
interaction (e.g., Too mre & Toomre|1972|[D'Onghia et al.|2010| >, 
and the ability to classify galaxies, either visually or using an au- 
tomated method, is strongly dependent on the surface brightness of 
the tidal features relative to the limitations of the data. With large 
data sets, such as those obtained from large surveys, the imaging 
is often not very deep, and the samples are large enough that vi- 



sual classifications quickly become impractical (although see Darg 
et al. 2010). As a result, many studies rely entirely upon the pro- 
jected separations as a metric for determining the merger phase of a 
galaxy pair, by correlating large r p with a large time elapsed since 
pericentric passage, and small separations with galaxies actively 
undergoing a close passage. The SFR enhancement trends with 
r p are thus usually interpreted as the signature of a sharp increase 
in SFR at pericentric passage, with the strength of the enhancement 
dwindling as the separations between the galaxies increases (e.g., 
Barton et al.|2000|[Barton Gillespie et al.|2003||Nikolic et al.|2004[ 



Woods et al.j2010p . This physical model would indicate that galax- 



ies are able to rapidly funnel gas to the central regions of a galaxy 
after an interaction, and that gaseous concentration is rapidly con- 
verted to stars, with the SFR declining as the galaxies separate and 
the burst ceases. 

To date, there has been no simultaneous statistical study of 
both the gas phase metallicity and the star formation rates in a sam- 
ple of galaxy pairs. Both the SFR and the gas phase metallicity 
are expected to change significantly throughout the merger as gas 
flows are induced through the galaxy; studying both of these quan- 
tities simultaneously allows us to gain insight into when changes 
are induced. In this work, we will tackle this issue by using a sam- 
ple of pair galaxies from the Sloan Digital Sky Survey Data Release 
7 (SDSS DR7; |Abazajian et aL|200"9") . Our sample will be used to 
precisely quantify differences between the pair and the control sam- 
ples over a wide range of projected separations and mass ratios. A 
suite of simulations based on the work in T orrey et al.| (2012 ) is de- 
veloped to allow us to analyse the trends with r p for characteristic 
markers of the merger stages. 

To this end, we use the sample of galaxy pairs in the SDSS 
DR7 compiled by Patton et al. 1(2011) and the metallicities calcu- 
lated in Scudder et al. (2012) to construct a clean pairs sample with 
stringent quality control measures. In Section|2]we review our sam- 
ple selection and define the control sample, along with a descrip- 
tion of our metallicity and SFR values. In Section [3] we describe 
our methodology for quantifying the changes in the SFR and the 
metallicity relative to the control sample and further quantify the 
changes for a subsample of morphologically disturbed galaxies. In 
Section|4]we describe the diagnostic statistics of the sample in more 
detail, and in Section [5] we compare our results to simulations in 
order to develop an interpretative framework for our trends with 
projected separation. We conclude with a comparison to previous 
work and an analysis of the physical picture our results suggest in 
Section|6] 

Throughout this work, we assume Sl M = 0.3, Q.a = 0.7 and 
H = 70 km s- 1 Mpc~'. 



2 SAMPLE SELECTION 

In this section we describe the criteria applied to the spectroscopic 
pool of galaxies in the SDSS DR7 from which we compile our final 
pairs and control samples. We wish to select a tightly controlled set 
of pair galaxies with a high probability of being a clean selection 
of physically associated and interacting systems. We require con- 
sistently calculated quantities across the pairs and the controls in 
order to be able to accurately compare between the samples. 

In order to qualify for the pairs sample or the control sample, 
a galaxy must have a stellar mas^Jand sufficiently strong emission 



' Unless otherwise stated, all mass values in this work are stellar masses. 
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lines so that a gas-phase metallicity and star formation rate can be 
calculated. In order to prevent the issues with Sloan photometry de- 
tailed in Simard et al. (201 1 1 from affecting the mass estimates, the 
masses used here have been recalculated from the MPA catalogu^] 
using their inferred mass to light ratios (M/L). The relationship be- 
tween the MPA-derived M/L and (g - r) colour is determined for 
the total MPA catalogue. This relationship is then used to derive 
masses from the g — r colours derived from the updated photom- 
etry presented in Sim ard et al.| ( |20lT|. With a n ew M/L value for 
each galaxy, the luminosities from Simard et al.| ( [201 lj > are used to 
calculate a mass without needing to rely on the Sloan magnitudes. 
We use the SFR values presented in the MPA catalogue, calculated 
according to the Brinchmann et al. ( 2004]) work. For galaxies with 
sufficiently high emission line signal to noise (S/N), Brinchmann 
|et al.| (2004| > uses a set of 6 emission line fits to the SDSS spec- 
tra to calibrate the SFR; lower S/N galaxies are calibrated with the 
strength of the Ha line. Given the high S/N we require for emission 
lines in our sample, all SFR values will be from the multi-emission 
line fits. The median lcr error on any individual SFR value in our 
sample is ~ 0.09 dex. The SFR values taken from Sloan's 3 arc- 
second fibre can be corrected to total values using models fit to the 
photometry of the galaxy outside the fibre, taking into account the 
fraction of light not contained within the fibre. A full description of 
the aperture correction methodology and bias testing is present in 
the Brinchmann et al. (2004) work. However, for this work, we use 
the fibre values only. 

We further require that each galaxy have a consistently cal- 
culated metallicity, following the criteria of Scu dder et al.| ( [2"0T2"| >. 
In order to ensure that all emission lines needed for the metallic- 
ity calibration are within the spectral range of the SDSS, we im- 
pose a lower redshift limit of z > 0.02. |Kewley & Ellison (20081 
demonstrate the problems of comparing between metallicity cal- 
ibrations; a systematically calculated diagnostic is crucial. The 
metallicity calibrations of Scudder et al. (2012} use the adaptation 
of the Kewley & Dopita ( 2002} method presented in jKewley & El-| 
jUson} 2008) . Briefly, a S/N > 5 is required in Ha, H/3, [On]/l3727, 
[Om]i4959,/I5007, and [Nn]i6584. Additionally, the Balmer ratio 
(Ha/H/?) S/N is required to be > 5. Galaxies must also be classified 
as star forming on the diagnostic diagram of |Baldwin eTaIJ ( [l98TT >, 
using the Kauffmann et al. (20031 diagnostic line. Using duplicate 
spectra, the standard deviation of the metallicity can be calculated, 
and results in a median error of 0.015 dex. We refer the reader to 
|Scudder et al.|p012[ ) for a complete description of the emission 
line quality control and metallicity calibrations. 

Our galaxy pairs sample is based upon the spectroscopic cat- 
alogue of |Patton et al.|p01 1) . For a full description of the sample 
selection algorithm, we refer the reader to that work. Briefly, |Patton| 
|et al.|p011) select galaxies with a spectroscopic close companion 
within r p < 80 kpc, with a velocity difference Av < 10,000 
km s~' , and within a mass ratio of 10: 1 . This results in a pairs sam- 
ple of 23,397 galaxies in pairs or higher order multiple^] Our mas- 
ter pairs catalogue is the result of a re-running of the algorithm of 
|Patton et al. (2011) with a minor alteration to the culling process 
meant to account for fibre collisions in Sloan. Due to the fixed fi- 
bre collision constraint of 55" ( |Strauss et al.|2002] l, galaxies with a 
projected companion nearer than 55" will be preferentially missed 



2 Available here: http://www.mpa-garching.mpg.de/SDSS/DR7/ 

3 ~8% of these galaxies have more than one close companion; 1.6% have 
more than 2 companions. Excluding these from our sample does not signif- 
icantly impact our results. 



in the spectroscopic sample. Some galaxies with separations less 
than the collision limit do have spectroscopic companions, primar- 
ily as a result of overlap between plates. In order to compensate 
for this incompleteness effect, 67.5% of galaxies in pairs with sep- 
arations > 55" are randomly culled from our sample, according to 
the incompleteness calculations of Patton & Atfield (2008 ). In Pat- 
|ton et ak] l |2011) , this cull was done without taking into account 
the culled galaxy's companion; i.e., a galaxy could remain in the 
sample while its companion was culled. For the sample used in this 
work, the culling has been repeated, but once one galaxy in a pair 
has been excluded from the sample, its companion is also removed. 

Galaxies with Av values at the high end of the distribution 
are unlikely to be physically associated. To select a sample of sys- 
tems which are more likely to be physically associated, a Av < 300 
km s~' limit is imposed to minimize the influence of projected pairs 
(Patton et al. 2000). This cut is in line with Av cuts in other works, 



which usually range between 250 - 500 km s '(e.g., |Lambas et al. 
[20031 pilson et al.||2008b[ |Patton et aLl[20lT) . 1,899 galaxies in 
the master pairs catalogue pass all criteria, and comprise our final 
galaxy pairs sampl^] 



2.1 Matching to Controls 

In order to determine whether an interaction leads to changes in 
galactic properties, we must also define a control sample. To this 
end, we define a sample of non-pair galaxies as any galaxy which 
meets the mass, SFR and metallicity requirements of the pairs sam- 
ple, but is not part of the full pairs catalogue (the "control pool"). 

The largest potential sources of bias between a sample of 
galaxies and its control are differences in the stellar masses, red- 
shifts, and environments of the galaxies in the two samples jPerez] 
|et al.|2009) . In order to eliminate these sources of bias from our 
sample, we wish to match our control sample in all three parame- 
ters, such that the normalized distribution of the final control sam- 
ple matches that of the pair galaxies. Simply matching in stellar 
mass and redshift could bias the results, as it has been shown that 
pairs tend to exist in higher density environments than non-pairs 
(e.g., |Barton et aT]|2007[ |Patton et aL]|2011), and there is a well 
known trend of SFR with density (e.g., Kauffmann et al.|[2004] 
|Poggianti et al.||2008| |Cooper et al.||2008) . Our density metric is 
calculated as an overdensity by taking the 5th nearest neighbour 
density and normalizing by the median density in a redshift slice of 
±0.01 around the galaxy. 

Galaxies in our pairs sample are therefore matched to galaxies 
in the control pool simultaneously in mass, redshift, and local den- 
sity, in a similar way to Ell ison et al.| p010). Briefly, our algorithm 
finds the galaxy from the control pool which is the simultaneous 
best match in all three parameters for each pair galaxy. After every 
pair galaxy has been matched to a control, a Kolmogorov-Smirnov 
(KS) test is used on the total distributions of pairs and control. If the 
KS test finds that the distributions are consistent with being drawn 
from the same parent distribution at > 30%, then the tentative con- 
trols become part of the control sample. Matching continues with- 
out replacement until 10 control galaxies have been matched to 
each pair galaxy in our sample, or the KS test returns a probability 
< 30%. For our sample of galaxy pairs, the algorithm ran without 
failing the KS criterion; the final control sample therefore contains 
18,990 galaxies. 



We do not require that both a galaxy and its companion pass all criteria. 
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Figure 1. Normalized distribution of the pair galaxies (solid blue) and the 
control sample (dashed black), matched in mass (top), redshift (middle), and 
density (bottom). Density is measured via the distance to the 5th nearest 
neighbour. Each pair has 10 controls. KS tests result in a probability of 
~ 99.8% that both the stellar mass and density distributions were drawn 
from the same parent distribution, and a probability of ~ 83.8% that the 
redshift distributions were drawn from the same parent population. 



3 SFR AND METALLICITY OFFSETS 

With a well-matched control sample of non-pair galaxies, it is now 
possible to proceed to compare the SFRs and metallicities of the 
pairs and controls. As we are using the fibre values for the SFRs, 
both our SFRs and metallicities are values for the central 3". Our 
adopted methodology is the same as Scudder et al. (2012 1, which 
itself is a modification of the method presented in Patton et al.| 
< |20 1 1| >. Each galaxy in a pair is matched to 10 controls, which all 
have metallicity and SFR values. The median SFR or metallicity 
of the 10 controls is taken, and subtracted from the value for the 
matched pair galaxy, as expressed in Equations [T] & [2] We define 
the resultant value, Alog(SFR) or Alog(0/H), to be our 'offset' val- 
ues from the control. 

Alog(0/H) = (log(0/H) + 12)^-^/2 (log(0/H) + l2) conlrols , (1) 

Alog(SFR) = log(SFR) ;)fl „. - A , 1/2 log(SFR) l . „„. , J , (2) 

where jii/2 signifies the median. 

Positive offset values indicate enhancements over the control 
sample, i.e., metallicity enrichment or SFR enhancement, whereas 
negative offset values indicate suppression relative to the control 
sample, i.e., metallicity dilution or suppressed SFR. This method 
has the advantage of allowing the quantification of the changes in 
SFR and metallicity on a galaxy by galaxy basis, relative to the 10 
control galaxies to which it is matched. Both SFR and metallicity 
are known to have strong mass dependen ces (jTremonti et al.|2004| 
|EUison et al.|2008a||Noeske et al.|2007l|Elbaz et al.|2011| >, but as 
pair galaxies and their controls are tightly matched in mass, a com- 
parison to the matched control galaxies is effectively a comparison 
at fixed mass. 

The calculated SFR and metallicity offsets are plotted 
against projected separation (r p ) in Figure [2] The distribution of 
Alog(SFR) is already visibly offset from zero. This trend is made 
clearer in Figure [3] which shows the same sample of galaxies, but 
binned in r p . SFRs are significantly offset from the control sam- 
ple out to 80 hzl kpc, the widest separation probed in our sample, 
with an increase in the offset magnitude at the smallest separations 
(r p < 15 fc^J kpc) to ~ 0.25 dex (about a factor of two). At separa- 
tions > 30 kpc, the offsets maintain a roughly constant magni- 
tude of ~0. 1 1 dex (a 30% enhancement). 

Alog(0/H) displays a much weaker trend with r p than the 
SFR, but similarly shows a sharp increase in offset magnitude at 
the smallest separations. Metallicity values are significantly sup- 
pressed by -0.02 dex out to ~ 60 tCj\ kpc. The plateau seen in the 
SFR offsets is not apparent in the metallicity offsets for the total 
sample of pairs. The relative weakness of this trend is not surpris- 
ing, as metallicity shifts due to interactions are usually found to be 
a small magnitude effect, generally of order -0.03 to -0.05 dex 
(e.g.,|Michel-Dansac et al.|2008||Cooper et al.|2008||EUison et al.| 
|2008b||2009l i~ 



The normalized distributions of stellar mass, redshift, and den- 
sity for the pairs and control samples are displayed in Figure[T| KS 
tests on the final distributions result in probabilities of 99.77% for 
the stellar masses of the two samples, 83.84% for the redshifts, and 
99.85% for the density distribution, indicating that the galaxy pair 
and control samples are very unlikely to have been drawn from 
different parent populations. Pairs and their controls are typically 
matched to within 0.001 dex in total stellar mass, 6 x 10~ 6 in red- 
shift, and 0.0008 dex in density, with a median range within the set 
of control galaxies of 0.07 dex in mass, 0.005 in redshift, and 0.09 
dex in density. 



3.1 Major & Minor Mergers 

It is interesting to consider the difference in SFR trends for major 
versus minor mergers. Our sample is large enough that we can di- 
vide according to the mass ratios of the interacting galaxies. We 
note that we divide our galaxies by stellar mass instead of by mag- 
nitude (cf. |Lambas et al.|2003|[Woods & Geller|2007||Lambas et al.| 
|2012| l; see |EUison et a l. (2008b I for a more extensive discussion 
of the use of luminosity, rather than mass, ratios. Galaxies within 
a mass ratio 0.33 < M host IM cm „ panjon < 3.0 are considered ma- 
jor mergers, M lmst IM companion > 3.0 are more massive companions 
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Figure 2. Offset values for all 1 899 galaxies in pairs sample as a function of 
projected separation. The top panel shows the SFR offsets, and the bottom 
panel shows metallicity offsets. The horizontal black dotted lines indicate 
the zero line. 
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Figure 3. All 1 899 galaxies in pairs sample. The top panel shows the SFR 
offsets, bottom panel shows metallicity offsets. Points are median values for 
the bin, and error bars are standard error on the median. Horizontal black 
dotted lines indicate the zero line. 



in a minor merger, and M hos ,IM companion < 0.33 are less massive 
companions in a minor merger. Our sample contains 1116 galax- 
ies in major mergers, and 783 in a minor merger, of which 184 are 
the more massive companion, and 599 are the less massive com- 
panions. The smaller number of more massive companions in mi- 
nor mergers is likely due to the increased probability of high mass 
galaxies hosting an AGN (e.g., Kauff rnann et al.|20 03 1 and thereby 
being removed from our star-forming sample. 

The trends in Alog(SFR) in major, more massive companion, 
and less massive companion subsamples with r p are plotted in the 
top, middle, and bottom panels of Figure|4]respectively. Grey back- 
ground points are the SFR offsets from Figure|3]in all panels. No- 
tably, galaxy pairs of all mass ratios are enhanced relative to the 
control at all separations. The shape of the offset vs. r p relation 
seen in Figure[3]is apparently driven by the major mergers. This is 
unsurprising, as major mergers are generally observed to show the 
strongest effects (e.g., | Woods et al.|2006| |Woods & Geller|2007"l 
Ellison et al. 2008b I, and they make up the majority (59%) of our 
sample. 

With the exception of the more massive companion bin at 50 
kpc, both the more massive and less massive companions in a 
minor merger show relatively Hat enhancements at the same mag- 
nitude (~ 25% enhancement) at wider separations. However, the 
innermost bin shows a different response between more massive 
and less massive companions in a minor merger. The less massive 
companion shows no strong enhancement at small separations, dis- 
tinguishing it from both the more massive companions and the ma- 
jor mergers, both of which show an increase in the magnitude of the 



SFR enhancement. More massive companions, on the other hand, 
show similar enhancement as the major mergers at small separa- 
tions, at 2 times stronger than the control; major pairs are enhanced 
by a factor of 1.9. However, as the more massive companions have 
the smallest sample, each binned point in Figure [4] only has 20-30 
galaxies. While the use of a median means that we are not biased 
by one or two outlying points, the poor number statistics may re- 
sult in an anomalously high point due to poor sampling of the total 
distribution. A more detailed discussion of the offset distributions 
is presented in j |4. 1| 

As the significance of the metallicity offsets is much weaker 
than that of the SFR offsets, splitting the metallicity offsets into 
mass ratio bins does not provide any additional information. 

3.2 Visual Classifications 

If the SFR enhancements and metallicity dilutions are truly be- 
ing driven by the tidal interactions of galaxies in pairs, then se- 
lecting a subsample which shows morphological evidence of a re- 
cent tidal interaction ought to amplify the effects seen in Figure [3] 
(e.g., |Michel-Dansac et al.|2008[|Lambas et al.|2012| which found 
stronger effects in a morphologically disturbed subsample). To this 
end, we select only those galaxies which show strong tidal arms 
or other asymmetries induced by an interaction. Although the pairs 
sample is designed to minimize the inclusion of projected pairs, 
physically bound galaxy pairs which have not undergone their first 
pass should have SFRs and metallicities close to the control, and 
will weaken the interaction-triggered signal. Any remaining pro- 
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Figure 4. SFR offsets, separated by mass ratio. Grey background points in 
all panels show the trends from the total sample (Figure[3j. The top panel 
shows the trend for all galaxies in major mergers (1116 galaxies). Middle 
panel shows the more massive companion in a minor merger ( 1 84 galaxies) 
and the bottom panel shows the trend for the less massive companion in a 
minor merger (599 galaxies). The more massive galaxies in a minor merger 
show an extremely strong enhancement at small separations, whereas the 
less massive companions show a much more consistent effect without any 
small r p increase in offset. Both the more massive and less massive com- 
panions show similar levels of offset at wide separations (~ 0.1 dex). 



jected pairs with low Av values would not be excluded from the 
sample, and, as physically dissociated systems, would also weaken 
the signal from the interacting systems. 

In order to eliminate the weakening effect due to either pre- 
interaction or projected pairs, all 1899 galaxies in the sample were 
visually classified for signs of morphological disturbances, as nei- 
ther sample would be expected to show morphological signs of 
recent interactions. Galaxies were flagged as either 'visibly dis- 
turbed' or 'not visibly disturbed'. Galaxies which fall into the cat- 
egory of not being visibly disturbed will include pre-interaction 
galaxies and interacting galaxies whose tidal features are below the 
surface brightness limit of the SDSS imaging, along with galaxies 
which are not truly interacting, so this category is not a useful diag- 
nostic on its own. However, the set of galaxies which falls into the 
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Figure 6. SFR and metallicity offsets for 1 105 pair galaxies flagged as dis- 
turbed. Grey background points are the trends from the total sample (i.e., 
Figure [3J. SFR offsets are of approximately equal magnitude in the dis- 
turbed and total samples, but the metallicity offsets are significantly larger 
in magnitude for the disturbed sample than for the total sample. The median 
value of Alog(0/H) triples in magnitude when the disturbed sample only is 
taken. 



disturbed galaxy classification should be a clean sample of galaxies 
which have already had a close encounter with a companion. 1 105 
galaxies are classified as disturbed, and the remaining 794 as 'not 
visibly disturbed' (See Figure[5]for 5 randomly selected examples 
of each classification). 

The SFR and metallicity trends with r p for the subsample of 
visibly disturbed galaxies are shown in Figure [6] Grey points are 
the total sample points from Figure [3] the disturbed subsample is 
overplotted in coloured points. The SFRs do not show a systematic 
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boost to their offsets in the morphologically disturbed sample; most 
points in the disturbed subsample are consistent with, or slightly be- 
low, the total SFR offsets. A KS test reveals an 89.78% chance that 
the disturbed and total samples were drawn from the same parent 
population. However, the median metallicity has shifted to more 
metal-poor values in all r p bins, in some cases by nearly -0.02 dex. 
The median Alog(0/H) has tripled in magnitude to -0.03 dex for 
the disturbed sample. The disturbed and total metallicity samples 
have a KS probability of being drawn from the same distribution of 
only 1.29%. Within this subsample, galaxies are significantly metal 
poor out to 80 kpc, the same range over which the SFRs show 
enhancements. 

The fact that the SFRs show no significant enhancement in 
the disturbed subsample is a puzzling contrast to the significantly 
lower metallicities. To ensure that the shifts between the total sam- 
ple and the disturbed subsample are not simply due to the smaller 
size of the disturbed subsample, we bootstrap a random sample of 
1 105 galaxies from the full pairs sample, and find the median value 
of the random subsample in both SFR and metallicity. This resam- 
pling was then repeated 20,000 times. We find that it is very easy 
to obtain a SFR offset of +0.14 dex through random sampling of 
the Alog(SFR) distribution, but extremely difficult to reach median 
Alog(0/H) values of -0.03 dex at random. The bootstrap test in- 
dicates that the lower metallicity in the disturbed subsample is un- 
likely to be a statistical fluke due to a smaller sample size. However, 
the median SFR offset in the disturbed subsample is consistent with 
the total sample, and no significant difference in the SFR offsets of 
the disturbed sample is seen. 

One potential interpretation is that the timescales over which 
galaxies are visibly disturbed are more strongly correlated with the 
timescales in which metallicity dilution is the strongest, and only 
weakly correlated with strong star formation enhancements. There 
is some suggestion from simulations that tidal features are visible 
for both a shorter period of time than the SFR burst, and at earlier 
times than the SFR enhancement jLotz et a l. 20081. If this is the 
case, then a sample of galaxies selected to be morphologically dis- 
turbed would not necessarily be expected to identify the galaxies 
with the largest Alog(SFR). 



4 DISTRIBUTIONS OF SFR AND METALLICITY 
OFFSETS 

While the median offset values as a function of r p prove useful as 
a metric of the typical merger, it is interesting to explore the range 
of metallicity and SFR offsets. The distributions will provide, for 
example, insight into how often the most extreme offsets occur, as 
well as how frequently no statistical change is seen. Furthermore, 
since simulations suggest that the two galaxies involved in a major 
merger ought to show similarly enhanced SFRs (Montuo ri et al.| 
|2010|[Torrey et al.|2012) , we can search for observational evidence 
of simultaneous triggering on a merger-by-merger basis. 

4.1 Offset distributions: Pairs vs. Controls 

The pairs sample is first split into two bins of r p . We define all 
galaxies with r p < 30 ICjk kpc as 'close pairs', and any galaxy 
pairs with r p > 30 hzl kpc as 'wide pairs'. The distributions of 
Alog(SFR) and Alog(0/H) are shown in Figures|7]&[8]respectively, 
with blue solid lines for the close pairs sample and red dashed 
lines for the wide pairs. The control sample is also divided by tak- 
ing all control galaxies matched to close pairs as the r p < 30 



kpc control sample, and all galaxies matched to wide pairs as the 
r p > 30 kpc control sample. Control offsets were calculated in 
a similar way as the pair offsets. Each galaxy in the control sam- 
ple was compared to the set of other galaxies matched to the same 
pair galaxy. As control galaxies were matched in sets of ten, every 
control galaxy has 9 galaxies of similar mass, redshift, and local 
density, to which it can be compared. We therefore take the me- 
dian of the 9 other control galaxies, and find the difference between 
that median value and the control galaxy selected. This difference 
(analogous to the calculations in Equation [TJ produces a control 
offset for every control galaxy in both SFR and metallicity. As the 
control galaxies are tightly matched, the control offsets ought to be 
centred around zero offset in both SFR and metallicity, with no sig- 
nificant shift away from zero or between the close and wide control 
samples. Figure s|7j& [8] show the control samples as the solid grey 
(close controls) and dashed black (wide controls) lines. The two 
control samples trace each other extremely well, and are centred 
around zero. 

The top and middle panels of Figures |7]& [8] offer two ways 
of viewing the distribution of offsets for SFR and metallicity. The 
top panel shows the offsets in a discrete way, with the overall form 
of the offset distributions for the two subsamples of both pairs and 
controls, normalized to the sample size. The middle panel shows 
the cumulative distribution of the same samples, which illustrates 
the differences between the distributions more clearly. Figures|7]& 
[8] show the total pairs sample rather than the disturbed subsample 
in order to maintain the number statistics of our sample. While the 
full sample will contain a small fraction of galaxies which are not 
physically associated, and a larger fraction of galaxies which have 
not yet interacted, this should only increase the number of galaxies 
with offsets near zero, and will not affect the distribution of offsets 
further from zero. 

In the top panel of Figure |7] it is already clear that both the 
close and wide pair samples are shifted systematically to higher 
offsets than their respective control samples, a trend which is even 
more dramatic in the middle panel. The cumulative distribution is 
calculated as the fraction of the total sample which has an offset 
greater than a given threshold. As the threshold values increase, 
the fraction of the control samples which have strong offsets de- 
creases; the pairs decline more slowly. KS tests confirm this vi- 
sual offset; both close and wide pairs are inconsistent with being 
drawn from the same parent population as their controls at > 6cr. 
To quantify this in a different way, the median Alog(SFR) of the 
close pairs sample is +0.21 dex (or a 60% enhancement over the 
control; see also Figure [3j, whereas the wide pairs sample has a 
median Alog(SFR) of +0.11 dex (a 30% enhancement). 

The bottom panels of Figure|7]shows the ratio of the pair frac- 
tion greater than a given offset (seen in the centre panel of Figure 
|7| relative to its respective control fraction, as a function of offset 
value. The black dotted line indicates a ratio of 1:1, i.e., that the 
pair and control samples have exactly the same fraction of the total 
sample at that offset or greater. This panel illustrates the relative 
frequency with which a given offset will appear in the pairs sample 
instead of the control. Values higher than one indicate that the off- 
set is preferentially found in the pairs sample. The shaded regions 
indicate the range of lcr errors, calculated from y/N statistics. This 
panel indicates that as Alog(SFR) becomes more extreme, it is in- 
creasingly likely to observe these offsets in the close pairs sample, 
rather than the controls or the wide pairs sample. The wide pairs 
sample also shows an excess of positive Alog(SFR) relative to the 
control sample, although only up to offsets less than +0.85 dex (7 
times the control value) Therefore, while the median SFR offset 
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is only an increase of 40%, there is an excess of close pairs with 
SFR enhancements up to a factor of 10 stronger than the control. 
Although only ~ 3% of pairs have excesses at this level, this is 3 
times more than exist in the control. Wide pairs, by contrast, only 
show excesses up to a factor of 7 (0.85 dex). 

Figure [8] is set up in the same way as Figure [7] The system- 
atic shift of the metallicities in the pairs sample is less pronounced 
than that of the SFR offsets, but is still visible. KS tests also verify 
this visual trend, with the metallicities for both close and wide pairs 
being inconsistent with being drawn from their parent distribution 
at > 4cr. Here the cumulative distribution is calculated as the frac- 
tion of galaxies with offsets less than the tested offset. The close 
pairs show a distinct shift to lower values relative to the control, 
whereas the wide pairs do not show a particularly strong shift. The 
difference in median is also much weaker; close pairs have a me- 
dian offset of -0.034 dex (8%), and the median for the wide pairs 
is -0.01 dex (2%). 

This distinction between close and wide pairs in the metal- 
licities is confirmed in the bottom panel, which shows the relative 
frequency of a given offset between the pair and control samples. 
Here we see that the the wide pairs do not show a significant excess 
of metallicity offsets relative to the control at any magnitude (see 
also Figure [3j. However, the close pairs show a significant excess 
of metal-poor galaxies at all negative offsets up to -0.25 dex (up 
to 1.78 times lower than the control). Offsets more extreme than 
-0.25 dex are so rare that the significance is overwhelmed by poor 
number statistics. The excess of very low metallicity galaxies is in- 
teresting, considering that the median offset for the pairs is 0.02 

dex, or 5% lower than the control. The discrepancy between the 
magnitudes of the excess offsets in the pairs sample and the overall 
median seems to indicate that the median offset is at least partially 
driven by a small number of galaxies which are very strongly offset 
from the controls. 

If, instead of splitting the total pairs sample into close pairs 
and wide pairs, we divide the sample into bins of mass ratio, we 
can investigate how the offset distributions depend on the mass ra- 
tio of the merger. The pairs sample is divided into major mergers 
(0.33 < yij] 0SI /Mcompanion ^ 3), less massive galaxies in a minor 
merger (M Ho „/Mc „,, Mn ,„„ < 0.33), and the more massive galaxies 
in a minor merger (M Hos ,/Mc om p an j,, n > 3). The control samples are 
split so that the control galaxies matched to pair galaxies in each 
bin in mass ratio are assigned to their respective bins. Figure [5] is 
constructed identically to Figures [7] & [8] The top panel of Figure 
[9] shows that all three mass ratio bins (in colour) are systemati- 
cally shifted to higher SFR values than their control samples (in 
grey). The control samples overlay each other reasonably well and 
are centred around zero, as expected. Given the similarity of the 
median values in Figure [4] it is unsurprising to see that the 3 dis- 
tributions are shifted by approximately the same amount relative 
to the control samples. KS tests give probabilities indicating that 
none of the galaxy pair subsamples are likely to be drawn from 
the same parent distribution as their respective control samples at 
> 5<r in all cases. The middle panel of Figure [9] shows the frac- 
tion of each mass ratio subsample which has an offset greater than 
a given threshold, as a function of that threshold. Here again we 
can see that the three distributions are roughly consistent in their 
shift to higher Alog(SFR), relative to the control sample. KS-tests 
on the distributions of pair offsets in the 3 mass ratio subsets indi- 
cate that the 3 distributions are consistent with being drawn from 
the same parent population; none of them are inconsistent with the 
null hypothesis at > 3cr confidence. 

The similarity between the offset distributions between mass 
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Figure 7. Top panel: histogram of SFR offsets for close pairs (blue solid), 
wide pairs (red dashed), and their controls. Control offsets are calculated 
between all controls matched to the same galaxy. Both close and wide pair 
galaxies are visibly shifted to higher offset values than the controls. Middle 
panel: cumulative distribution of SFR offsets, indicating the fraction of the 
galaxies in the sample with offsets greater than a given value. Galaxies are 
divided in the same way as the top panel. At all offset values, the close pairs 
sample shows a higher fraction of galaxies with strong SFR enhancement. 
The bottom panel shows the ratio of the fraction of pair galaxies above 
the threshold offset to the fraction of the control above that offset, with the 
shaded region indicating ^fN errors. Black dotted line shows a ratio of 1. 



ratio bins is reinforced in the bottom panel of Figure [9] All three 
samples show an excess of positive values of Alog(SFR) relative to 
the control. Both major and minor mergers seem equally effective 
at inducing offsets up to 0.45 dex above the control value, or an en- 
hancement of a factor of ~ 3. These intermediate offsets are roughly 
1.3 times as likely to occur in the pairs sample as the control in all 
samples, with no significant distinction between the more massive 
and less massive companions in the minor mergers, or between the 
minor mergers and the major mergers. SFR enhancements of a fac- 
tor of 2 (+0.3 dex) occur in 30-35% of the pairs sample, and are 
1.35-1.65 times more likely to occur in the pairs than in the con- 
trol. However, offsets more extreme than 0.8 dex are reached al- 
most exclusively through major mergers. These offsets are just as 
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Figure 8. Same as FigurejT] but for metallicities. The top panel shows the 
distribution of the pairs and control samples, split into r p < 30 kpc and r p > 
30 kpc bins, coloured same as Figure^] Middle panel shows the cumulative 
distribution of metallicity offsets, and bottom panel shows the ratio of the 
pair and control fractions for > 30 and < 30 kpc bins. The shaded region 
indicates ^/N errors. 

rare as before, occurring in roughly 5% of galaxies, but are more 
than twice as likely to occur in major pairs as in the control. 

In summary, extreme offsets are rare in the pairs sample, but 
are most likely to be found in the close pairs with approximately 
equal masses. SFR enhancements of + 1.0 dex (a factor of 10 higher 
than the control) are 3 times as likely to occur in the pairs as in the 
control, but occur in only 3% of the total pairs sample. Similarly, 
metallicity offsets of -0.25 dex (a factor of 1.78 lower than the con- 
trol) are 1.5 times as likely to occur in the pairs as in the control, but 
also are only found in 3% of the pairs sample. Conversely, modest 
SFR enhancements, up to a factor of ~ 3 beyond the control, can be 
found as readily in the minor mergers as the equal mass pairings. 

4.2 Evidence for synchronised SFR triggering 

Previous work has found that galaxies in pairs often have both 
galaxies show enhanced star formation more frequently than would 
be expected at random ( |Xu et al.|2010) . This is interpreted as the 
signature of the two galaxies in an interaction undergoing syn- 



Figure 9. Same as FigurejT] but split by mass ratio instead of r p . The top 
panel shows the distribution of offsets for the three mass ratio bins: major 
mergers (0.33 < MHastft^Companion ^ 3), the less massive companion in 
a minor merger (MHostP^Companian < 0.33), and the more massive galaxy 
in a minor merger (M/y os/ /Mc ,„ /M ,„ „ > 3). The control galaxies for the 
galaxies in each mass ratio bin are plotted in grey. The centre panel shows 
the fraction of each of the three samples (and respective controls in grey) 
which has an offset greater than a given value. The bottom panel shows the 
ratio fraction of the pairs galaxies above the threshold offset to the fraction 
of the control above that offset. Shaded regions indicate errors. 



chronised SFR enhancement, as both galaxies undergo similar tidal 
torques due to their companion. This physical picture is supported 
by simulations of equal mass mergers, where both galaxies tend to 
show similar responses to the interaction (e.g., |Torrey et al.|2012| >. 
Some observational evidence for merger-driven synchronicity be- 
tween galaxies has already been found in galaxy pairs with AGN. 
|EUison et al.H20lT) found that in galaxy pairs at small separations, 
a galaxy hosting an AGN was twice as likely to have a companion 
also hosting an AGN than would be expected at random. We there- 
fore wish to investigate whether the Alog(SFR) values also show 
evidence of synchronised enhancements. 

We first select only those galaxies from the total pairs sam- 
ple where both the pair and the companion are found in our sam- 
ple. This reduces us to 45% of our original sample; 425 galaxy 
pairs remain in the sample (850 galaxies). We then plot the host 
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will have the same distribution of offsets as our true pairs sample. 
We find that the scrambled pairs populate the double enhancement 
quadrant with 40.99% of the total sample. 13.28% of the scrambled 
pairs sample falls in the double-suppression quadrant, identical to 
the true pairs sample. However, the double enhancement quadrant 
has 17% more galaxies in the true galaxy pair sample than in the 
scrambled pairs. 

The SFR offsets have some intrinsic scatter, and galaxies very 
close to the zero line may simply be scattered to one side of zero 
or another. To minimize the effect of this scatter, we can include 
a buffer around the zero value, such that galaxies within a certain 
range around the offset value are not counted in this fractional 
counting. If the buffer is used, the strength of the double SFR en- 
hancement increases. If galaxies which have offsets beyond ±0.13 
dex (the 25th percentile for the control distribution) are considered 
in the fractional calculation, then the true pairs are 34% more likely 
to show double enhancements. 



Figure 10. For the subset of galaxy pairs where both the host and the com- 
panion are in our sample, (~45% of the sample; 425 galaxy pairs), the dis- 
tribution of the galaxy pair Alog(SFR) vs. its companion's Alog(SFR). The 
colour bar indicates the projected separations of the pairs. 48% of the sam- 
ple is found where both the pair and its companion have Alog(SFR) > 
(top right quadrant). Using a random galaxy pair match of control sample 
galaxies, we find each quadrant to be equally populated. Galaxies in pairs 
are 1.9 times as likely to show double SFR enhancements than a random 
control sample, and 0.6 times as likely to show double suppressions (bot- 
tom left quadrant). Relative to a false pairs sample, the galaxies in true pairs 
are >17% more likely to show double enhancements. 

Alog(SFR) versus its companion's Alog(SFR) (Figure [lOf, where 
the points are colour-coded according to projected separation. Each 
galaxy pair is plotted only once on this diagram. It is clear that 
there is an overabundance of points in the double-enhancement (top 
right) quadrant of the figure; this quadrant contains 48.0% of the to- 
tal pairs sample. The double-deficit quadrant contains only 14.59% 
of the sample. To gain a sense of what would be expected at ran- 
dom, we pair random control galaxies with each other, and find 
the distribution of points for the control offsets calculated in Sec- 
tion|4] As expected, each quadrant contains 25% of the control pair 
distribution. This indicates that the pairs are 1.9 times as likely to 
fall in the double-enhanced quadrant than would be expected from 
the control. Splitting this sample into galaxies in major mergers 
and those in minor mergers results in an almost identical fraction 
of galaxies with doubly enhanced SFRs (47.42% in major mergers 
vs 49.25% in minor mergers). The minor mergers show a slightly 
lower fraction of galaxies with double suppressions (2%) relative 
to the major mergers. 

However, since the Alog(SFR) values are higher in the pairs 
sample than they are in the control sample, simply having an ex- 
cess of galaxies with doubly enhanced SFRs does not necessarily 
indicate that there is synchronous triggering. This could instead be 
the result of comparing between two samples which are not drawn 
from the same parent distribution, rather than indicating that pairs 
of galaxies are showing correlated SFR enhancement. In order to 
eliminate the issue of the differing parent distributions, instead of 
comparing the population of the pairs in this diagram to the con- 
trol sample, we compare to the pairs sample itself. We reassign 
each galaxy in the sample to a random companion, thereby scram- 
bling the pairs sample. This results in a set of uncorrected pair 
Alog(SFR)s which can be plotted on the same diagram, but which 



5 COMPARISON WITH THEORETICAL MODELS 

The SFRs and metallicities of our pairs are significantly offset 
from the control over significantly larger distance scales than has 
been previously seen (e. g.,|Lambas et a l. 2003 ; Alons o~t al.|2006j 
|Woods & Geller|[2007) |Ellison et al.||2008b| >. The interpretation 
of the smoothly declining trend with increasing r p has tradition- 
ally been that galaxies promptly undergo a burst of star formation 
as the galaxies reach pericentre. This burst would then decline to 
a fiducial value as the galaxies separate. However, Figure [6] indi- 
cates that galaxies are still offset from the control values at ~80 fcjl 
kpc separations, and the wider separations show a plateau instead 
of a smooth decline to control values. The existing interpretation 
is insufficient to account for the wide separation plateau, as even 
a long lasting starburst should still show a smooth decline to the 
control values. We therefore turn to theoretical models to aid in the 
interpretation of these trends. 

We make use of simulations of major galaxy mergers devel- 
oped in |Torrey et al.| ( (20T2[ > to help interpret our results. A detailed 
analysis of varying mass ratios and orbital parameters has previ- 
ously been explored elsewhere (e.g., Di Matteo et al. 2007] |Cox| 
|et al.|20 08 1. Here, the goal is to take a general look at the SFR and 
metallicity changes in a merger as a function of r p . We refer the 
reader to |Torrey et al.| l (2012) for a detailed description of the sim- 
ulation setup and parameters. Briefly, Torrey et al. (2012) presents 
a set of N-body/Smooth Particle Hydrodynamics simulations us- 
ing Gadget-2 (Springel 2005). These models include cooling, star 
formation, feedback, and chemical enrichment. Within these simu- 
lations, the metallicities are defined as the mass-weighted average 
of the metallicities of all gas particles within a sphere of 1 kpc 
around the centre of the galaxy. Increasing the radius of this sphere 
only mildly alters the results of the simulations, and does not af- 
fect our conclusions. (As a comparison, the physical diameter of 
the SDSS fibre is generally of order of a few h^l kpc.) Galaxies are 
shown to be stable (i.e., do not develop a bar) for at least 2 Gyr 
when modelled in isolation. For the purposes of constructing an 
interpretive framework for the SDSS data presented here, we use 
a suite of 16 galaxy mergers, varying only the galaxy orientations 
between mergers. As a result, their orbits are kept constant through 
the suite. The model galaxies are merged with identical copies of 
themselves in all 16 simulations 



5 Our model galaxy is Disk B in Table 1 of|Torrey et al.| |2012|, but with 
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An example of the evolution of one of the 16 mergers is shown 
in Figure [TJJ This figure shows the time evolution of the sepa- 
ration of the two nuclei in the top panel, the change in metallic - 
ity in the centre panel, and the change in SFR for both galax- 
ies in the bottom panel, where t = is scaled to coalescence. 
Snapshots of the galaxies' gas densities are shown as insets at 
the top of the figure. Comparable values to the observational 
Alog(0/H) and Alog(SFR) values are extracted from the simula- 
tions. Alog(0/H) and Alog(SFR) are calculated as the difference 
between the merging galaxy's SFR or metallicity at any given time 
and the SFR or metallicity of a model quiescent disk at that same 
time. The two lines in the bottom two panels show the responses of 
the two galaxies in the merger. For major mergers, the close track- 
ing of the two lines is a general feature of the merger simulations. 
The magnitudes of the SFR enhancement and metallicity dilution 
may shift slightly from merger to merger, but the overall shape of 
these tracks with time is consistent across the suite of simulations. 
Therefore, in some runs, the metallicities may re-enrich past the 
initial value, but this re-enrichment is almost always bracketed by 
periods of metallicity dilution. Similarly, the strength of the initial 
SFR burst may vary slightly from merger to merger, but always di- 
minishes again, and is dwarfed in magnitude by the peak in SFR as 
the galaxies reach the end of their merger. 

Of particular note in Figure [TT] is the time delay between first 
passage and the strongest metallicity dilutions and highest SFR en- 
hancements after pericentric passage. Metallicity dilutions do not 
reach their strongest values until ~ 200 Myr after first passage; 
the peak in SFR enhancement is slightly longer, at ~ 275 MyiQ 
This delay means that the galaxies have had time to progress out 
to wider separations by the time the triggered SFR enhancement 
or metallicity dilution will be strongest. There also is a period of 
metallicity re-enrichment after the first starburst, prior to coales- 
cence. Recently formed stars will return enriched gas to the inter- 
stellar medium as they reach the end of their lifespans, which drives 
the galaxies to nearly return to their initial metallicity. The flow 
of metal-poor gas to the central regions of the galaxies continues 
throughout the interaction, re-diluting the central metallicities. As 
galaxies progress into final coalescence, metallicity offsets drop to 
their lowest values, and Alog(SFR) reaches the strongest enhance- 
ment of the merger sequence, due to the large torques exerted as 
the nuclei coalesce. 

In order to more directly compare the results of the simulation 
to our data, we must fold the Alog(0/H) and Alog(SFR) tracks into 
a plot as a function of separation, rather than as a function of time 
elapsed. Figures [T2fr and|13fr show metallicity and SFR tracks re- 
spectively for one galaxy's response to a merger, as a function of 
true physical separation. The metallicity and SFR are measured at 
fixed time steps (10 Myrs) throughout the merger, and points are 
colour-coded as a function of time (from dark blue at early times to 
yellow at late times). Figure [T"2"fr shows that the maximal dilution 
in Alog(0/H) after the first close passage seen in Figure [TT] occurs 
at a separation of ~ 60 kpc (in this particular merger), with the final 
metallicity drop due to coalescence seen at real separations of < 15 
kpc, as the galaxies coalesce. Similarly, the strongest SFR enhance- 
ment after first passage occurs in a relatively narrow peak between 
45-65 kpc, with the coalescent peak appearing at separations < 15 



a 25% gas fraction, and merged on all 16 orientations in Table 2 of |Torrey] 
[etaL](20T2) . 

These timescales are primarily set by the free-fall timescales of the sys- 
tems, given a rotational velocity of the simulated galaxy of 130 km/s. 



kpc. Interestingly, between the two peaks, the SFR appears to drop 
to nearly its original value; high SFR offsets are generally present 
at wide separations after a close passage, or at very small separa- 
tions due to coalescence. However, the narrowness of the peak in 
SFR at 60 kpc in separation space is a result of the coincidence 
of the galaxies being at apocentre during the period of time when 
the SFR enhancement strengthens. There is no reason to expect the 
apocentre and the SFR burst to coincide; this is a coincidence of 
our merger orbital parameters. If the galaxy takes longer to reach 
apocentre, the SFR enhancement would begin prior to apocentre, 
and this SFR peak would have a broader distribution in separation 
space. 

Physical separations are not an observable quantity, so in order 
to place the results of the simulations on more comparable footing 
with the data, we introduce the full suite of 16 mergers, and convert 
the physical separation tracks into projected separation space. The 
conversion takes each point from the simulations and multiplies by 
a random viewing angle in 3-dimensional space, defined as |cos(0)|, 
where cf> = sin (R) and R is a different random value between and 
1 for each point in the simulations. The cosine of a random angle 
accounts for the 2 dimensional spin of the viewing angle, and to 
account for the distribution of angles over the surface of a sphere, 
we use the arcsine of a random number between and 1, which 
should properly account for the 3 dimensional distribution of possi- 
ble viewing angles. We can then plot the simulated Alog(0/H) and 
Alog(SFR) as a function of projected separation. The contours of 
the scattered points are shown in Figures |12b & |13b respectively. 
The projected separation contours are time-weighted, as each point 
in the real separation simulations is measured at fixed time inter- 
vals. In regions of real space where the galaxies spend most of their 
time, there is a corresponding increase in the number of points at 
that distance. 

Strikingly, these contours outline the same general form as is 
seen in the SDSS sample of pair galaxies. The metallicity contours 
outline a low level plateau of metallicity dilution out to 70 kpc, 
with a sharper drop at very small separations. The contours do not 
show a marked change from the individual tracks, other than an 
expected blurring of minor features due to projection effects and 
the scatter introduced by including several merger orientations. The 
SFR enhancements undergo a more dramatic transformation be- 
tween real and projected space. The reasonably tight peak in SFR 
at wide separations has scattered to smaller separations, resulting 
in a distinct plateau of enhanced SFR out to 70 kpc. Converting 
to projected separations will scatter galaxies to smaller separations, 
so the wide separation peak has filled in the previous gap between 
the high r p post-pericentre peak and the peak due to coalescence. 

The cutoff of the offset plateau at ~ 70 ft™ kpc is a result of 
our simplified merger suite. As we do not explore a large range of 
potential orbits, varying only the galaxy inclination, none of our 
mergers reach separations beyond ~70 kpc after their first pass. 
With a wider suite of simulations, we expect this low level plateau 
would naturally extend as galaxies reach larger distances after their 
first pass (Patton et al„ in prep). Observationally, this would cor- 
respond to galaxies in either more weakly bound interactions than 
those we have simulated here, or galaxies with higher initial an- 
gular momenta. The residual peak at 60 ft™ kpc is also a function 
of the fixed initial orbital energy and angular momentum in all 16 
of our simulated mergers; keeping these two parameters the same 
results in fixed orbits. With a wider range of orbits, this feature 
should weaken into a flatter plateau, as galaxies with smaller and 
wider apocentres blend together. Since the strength of the tidal in- 
teraction is expected to directly impact the strength of the induced 
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Figure 11. Simulation snapshots, Alog(SFR) (within the central 1 kpc), & Alog(0/H) as a function of time. The dashed line marks first pericentre, and the dot- 
dashed line indicates coalescence. The top series of panels shows the surface density of the gas within the galaxy at different points in the merger. Alog(0/H) is 
the difference between the metallicity of a quiescent disk and the metallicity of the merging galaxy at any given time step, and Alog(SFR) is the difference 
between the current SFR at a given time and the SFR of a galaxy modelled in isolation. The two solid lines in the metallicity and SFR panels indicate the 
responses of the two galaxies in the simulation. The separation panel shows the intergalactic distance as a function of time, which remains fixed through the 
full suite of simulations; Alog(0/H) andAlog(SFR) change between simulations as the galaxy orientations change. 



SFR (e.g., |Di Matteo et al.|2007l|Cox et al.|2008| >, we also expect 
that with a broader range of tidal interactions, we should see a cor- 
responding increase in the range of SFR enhancements induced in 
the simulations. 

The qualitative similarity between the contour plots result- 
ing from this simple suite of 16 major mergers and the observa- 
tional results makes the framework of these simulations an appeal- 
ing one for the interpretation of the observed trends. In this sce- 
nario, the innermost peak in offset values should be due almost 
entirely to galaxies approaching coalescence, while the wide sepa- 
ration plateau is due to galaxies which have gone through a close 
passage, and are only showing the SFR and metallicity response 
at wider separations. These galaxies could ultimately merge, or be 
part of a population of fly-by encounters, which are also expected 
to show SFR enhancements (e.g., Di Matteo et al. 2007 ; Montuori 
|et al.|2010| >. 

It is necessary to caution that these simulations are not in- 
tended to function as a direct quantitative comparison to the data, 
but simply as a theoretical framework to help interpret the form of 
the signal observed in the SDSS data. With a basic set of simula- 
tions, we have not reproduced (or indeed, attempted to generate) a 



representative sample of the range of mergers that exist within the 
SDSS sample. The galaxies in the simulations are major mergers 
only, with a single progenitor galaxy, on a constant merging or- 
bit. We can state only that the simulations do not lie in a region 
of parameter space devoid of points in the SDSS data. With these 
concerns in mind, we do not overplot the SDSS data on the sim- 
ulation contours, nor do we encourage a direct comparison of the 
magnitude of these effects. There are a large number of parameters 
which could change the magnitude of the signal seen in our suite of 
simulations (e.g., a varying fraction of interloping or pre-pericentre 
galaxies, the influence of minor mergers, gas fractions, or initial 
orbits), many of which have been the subject of other simulation 
studies (e.g., Di Matte o et al.|2007||Cox et al.|2008) . However, the 
comparison of the general form of the trends with r p seen in the 
Alog(SFR) and Alog(0/H) figures is robust. The fact that the trends 
seen in the SDSS sample is reproducible using a straightforward 
set of simulations is encouraging. This match in form indicates that 
the simulation does not need to be tuned to the particulars of our 
sample to observe the same general trends. 
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Figure 12. Panel a: 1 sample track of the metallicity offset as a function of true separation. Colour indicates the progression of time through the merger, with 
dark blue at early times to yellow at late times. Each dot indicates a fixed time step of 10 Myrs as the simulation progresses. Panel b: 16 tracks, converted 
into projected separations. Colours indicate the density within the contours, with yellow indicating the highest density of points from the simulation and blue 
indicating the lowest density. The solid black line overlaid on the contours indicates the median value, and the dashed black lines indicate the 25''' and 75''' 
percentile range. 
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Figure 13. Same as Figure [12] but for SFR. The contours in Panel b outline the same general form as is seen observationally, with a wide separation plateau, 
and higher SFR enhancements at small separations. 



6 DISCUSSION 

We find that in a strictly selected sample of galaxy pairs from the 
SDSS DR7, the star formation rates are typically enhanced by at 
least 30% out to separations of 80 /z^j kpc. The metallicities are 
suppressed by -0.02 dex (~ 5%) within 60 kpc. We visually 
classify all galaxies in our sample for signs of morphological dis- 
turbance to further clean the sample of galaxies which have not yet 
undergone an interaction or of the remaining fraction of interloping 
galaxies. When only the disturbed subsample is taken, the metallic- 
ity trend increases in significance over the entire range of r p and is 
offset from the control values out to at least 80 fcj kpc by -0.03 



dex. Within the inner 30 kpc, the disturbed sample is offset by 
-0.04 dex (9%) in metallicity, and enhanced by 65% in SFR, rel- 
ative to the control sample. Although previous studies have found 
enhanced SFRs and diluted metallicities in samples of close pairs 
(e.g., |Ellison et al . 2008b, Kewley et al. 2006i), this is the first time 
that the changes in SFR and metallicity have been studied simulta- 
neously as a function of projected separation. 
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6.1 Star formation rate enhancements out to 80 kpc. 

Most previous work using large sample statistics has found that 
the enhancements seen in SFR are restricted to within 30-40 h^i 



kpc separations (e. g., |Xu & Sulentic|[T99"T| |Lambas et al.|[2003| 
|Alonso et al.|2006||Eilison et al.|2008b[|Robaina et al.|2009| l.lri- 
deed, in our sample, the strongest signals are present at the small- 
est separations. However, we find statistically significant offsets in 
metallicity and SFR out to the maximum separations of our sample 
(80 k£ kpc). 

Many previous works have used a mean control value as the 



point of comparison for their SFR enhancements (e.g., Lambas 



et al. 2003, Alonso et al. 2006; Ellison et al. 2008b; Lambas et al. 



2012[ >, whereas we have used "offsets" for individually matched 
galaxies. We repeat our analysis using the mean of the control 
sample as the baseline for our SFR enhancements. Our results are 
broadly unchanged; we still see significant SFR enhancement at 80 
fcl kpc. Evidently, the change in methodology is not entirely re- 
sponsible for the increase in sensitivity to changes in SFR at wide 
separations. If the methodology is not the dominant improvement, 
then the other major set of differences between our work and pre- 
vious works is in the definition of a control sample. While most 
of the previous works match their pair and control samples in red- 
shift, few of them are matched in stellar mass. As the SFR of a 
galaxy is strongly dependent on mass, if the distributions of mass 
between the pairs and the control galaxies are not well matched, 
this could well influence the strength of the overall SFR enhance- 
ment measured. For instance, |Lambas et al.| ( |2"0"0"3"] > and |Darg et al.| 
( 20 1 ) match to a control sample only in redshif t. Other works have 



used magnitudes as a proxy for mass ( |Woods & Geller|2007||Wong| 
[eT^[20TT] |Alonso et alT2006*l >. However, |Ellison et al.| (|2008b I 
showed that this method can dilute trends with galaxy properties 
with a mass dependence, such as SFR. The conversion from mag- 
nitude ratios to mass ratios is not 1:1 (see their Figure 3), so tak- 
ing magnitude ratios as a direct proxy for mass will introduce a 
significant fraction of galaxies of different masses than intended. 
|Perez etal.l (2009) compares different methods of selecting a con- 
trol sample, and finds that matching only on luminosity and redshift 
introduces significant biases into the result, resulting in artificially 
redder, more massive pairs, which live in higher density regions of 
space, relative to their controls. The biases introduced by match- 
ing in luminosity and redshift can be reduced by 70% by instead 
matching a control sample in mass, redshift, and local density. 

Our pairs sample is matched to a control in total stellar mass, 
redshift, and in local density. [Patton et aT] ( |201 1) and |Ellison et al.| 
(2008b I also match in mass and redshift. The extent of the SFR off- 
set trend is remarkably similar to the results found in Patton et al. 
(2011 1. The Patton et al. (2011 1 study used a superset of our sam- 
ple, as it did not require the gas phase metallicity calculations, and 
uses a comparable offset methodology to ours; direct comparisons 
between the results are robust. Patton et al. ( 201 1 ) find that the cen- 
tral (g-r) colours of blue cloud galaxies are consistently bluer than 
the control values out to 80 fcjl kpc (see their Figure 15). Ellison 
|et al.| ( [2008b| > sees a hint of this trend in major mergers, but the er- 
rors are large enough that the trend is not systematically significant 
to the widest separations, and the full sample only shows signif- 

kpc. Our sample is moderately 



icant enhancement at r„< 



30 h£ 



larger than that of the Ellison et al. (2008b) work, at 1899 galaxies 
versus their 1719, and we have a slightly more stringent Av require- 
ment of 300 km/s instead of 500 km/s, both of which would help 
in reducing error bars. Furthermore, Ellison et al. (2008b I use total 
star formation rates, rather than the fibre values used in this work. 



|Patton et al. ( 201 i| determined that the strong offsets in colour are 
primarily a nuclear effect, as the total (g - r) colours of the galaxy 
showed smaller offsets across all r p by a factor of ~ 4, and lacked 
the small separation peak to bluer colours. We use the aperture cor- 
rected SFRs provided in Brinchmann et al. (2004 Qto test whether 



we see an analogous weakening of the magnitude of the SFR off- 
sets. Consistent with the |Patton et al.| l(201 1 1 work, we find that the 
galaxy's total SFR shows a weaker trend than the fibre values, and 
lacks the central spike to higher SFR. However, SFR offsets for the 
total values are still visible at 80 hzl kpc, so while using the fibre 
values increases the magnitude of our SFR enhancements, the use 
of the fibre values alone is not enough to account for the new sen- 
sitivity. A similar comparison is impossible for the metallicities, as 
there is no correction from fibre metallicities to total metallicities. 
The nuclear concentration of the SFR enhancement indicates that 
it is much more likely to detect strong SFR offsets when looking 
at the fibre values rather than the total values, as in [Ellison et al.| 
l |2008b| >. We therefore propose that one of the main factors in in- 
creasing our sensitivity to small effects in SFR and metallicity is 
due to our tightly mass-matched control sample. Further aiding our 
sensitivity is the use of fibre SFRs rather than the aperture corrected 
values. 



6.2 Mass Ratios 

As suggested by many previous studies, the form of the trend with 
projected separation is primarily driven by the contribution of ma- 
jor mergers (e.g., [Woods et al.|2006| |Woods & Geller|2007| |HlT"| 
son et al. 2008b ). Investigating the more massive and less massive 
companions in minor mergers (mass ratios more extreme than the 
3:1 major merger criterion), we find that both the less massive and 
more massive companions show SFR enhancement at all r p . The 
more massive galaxies distinguish themselves from the less mas- 
sive galaxies only in the smallest separations (see Figure]?]). Major 
mergers have slightly higher median offsets, but the median value is 
very comparable to that of the minor merger enhancements. |Lan>] 
bas et al. ( 20 1 2 1 find a similar effect; minor mergers show enhanced 
SFRs at all masses, relative to the control. Their sample of major 
mergers is found to have SFRs enhanced at a slightly stronger level 
than the minor mergers, but in large part the offsets of the two sam- 
ples are consistent within error bars. We find that minor mergers are 
just as effective as major mergers at inducing SFR enhancements of 
less than a factor of 3 over the control. However, the rarest, most ex- 
treme starbursts occur almost entirely within major mergers. Since 
these extreme starbursts are so rare, our overall offset distributions 
are consistent between major mergers and minor mergers. 

Previous studies of SFR enhancement as a function of mass 
ratio have drawn conflicting interpretations from their data (e.g., 
|Lambas et a l. 2003 determined that the more massive galaxy in 
a minor merger is more strongly affected, while Wo ods & Geller| 
|2007| found the inverse), a closer inspection of the data indicates 
that these two previous observational results and that of our cur- 
rent work are not inconsistent. Woods & Geller ( 2007 1 found no 
evidence for a correlation between small r p and high SFR for the 
higher mass companion in a minor merger, whereas there was some 
evidence of correlation for the less massive companion, and con- 
clude that the less massive companion is more strongly affected by 



7 |Brinchmann et al"]{2004) use model fits to the photometry outside the 
fibre to make the correction from fibre to total values. A complete discussion 
of potential biases is presented in that work. 
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a merger than the high mass companion. However, since the SFRs 
of the galaxies in pairs are not directly compared to a control sam- 
ple, it is impossible to judge whether or not these galaxies are also 
systematically enhanced over the control value. For example, our 
sample of less massive galaxies in a minor merger shows system- 
atic SFR enhancement, but no correlation with the r p of the galaxy 
pair, so on this point our data do not necessarily conflict. 

|Woods & Geller| 12007) also found that when comparing the 
distributions of specific star formation rates (SSFRs), the distribu- 
tions of both the more massive and less massive galaxies in minor 
pairs were statistically consistent with the field sample. |Lambas| 
|et al. | (2003;) also find that galaxies in a minor merger show no sig- 
nificant SFR enhancement at any r p , compared with the average 
control SFR value. In contrast, we find that the distribution of star 
formation rates at a given mass are statistically different for galax- 
ies in minor pairs and our control sample; KS tests indicate that the 
less massive galaxies are unlikely to be drawn from the same parent 
distribution as the control at ~ 6cr, and at > 5c for the more mas- 
sive galaxies in a minor merger. We note that both Lamba s~et al.| 
l [2003) and |Woods & Geller| ( [2007) are subject to the same issues 
with matching a control sample in magnitudes instead of using stel- 
lar masses described in Ell ison et al.| ( |2008b") , which will weaken 
the sensitivity of their measurements. 

While our results are not in conflict with the observational 
results, our results do not appear to align with the expectations 
from simulations. There are relatively few simulations which have 
investigated the SFR enhancements of galaxies in unequal mass 
mergers. |Cox et al.| ( [2008l > find that high mass galaxies in unequal 
mass mergers are less likely to be tidally perturbed and drive strong 
SFR enhancements when interacting with a low mass companion, 
whereas the low mass companion will be strongly perturbed by a 
massive companion. Both |Cox et aL] (2008) and |Mihos & Hern-| 
|qu!st| ([T994J find that minor mergers can drive some gas inflows 
in the massive companion, but the inflowing gas is not necessarily 
converted into stars. Mihos & Hernq uTst| \ 1 994) find that the mas- 
sive companions only show enhanced SFR at coalescence, which 
does not help explain our trends at wide separations. |Cox et aT 



2008 1 is meant as an improvement upon the |Mihos & Hernquist 



1994 1 work, and finds that the massive galaxy is unlikely to un- 



dergo a starburst except in very specific cases. If this theoretical 
model were borne out, we might then expect to see significantly 
stronger Alog(SFR) in the less massive companions, particularly at 
small separations, where the galaxies are likely to be in the final 
stages of a merger. In contrast, we find that the more and less mas- 
sive galaxies show similar levels of SFR enhancement over most of 
the range in r p probed by our sample, with the more massive com- 
panions displaying higher SFR enhancement at the smallest sepa- 
rations. 



6.3 Magnitude of the SFR enhancement 

We can next compare the magnitude of our SFR enhancements to 
those found in previous works. We use a series of statistical tests 
beyond the median values to determine which magnitude of offsets 
are preferentially found in the pairs sample instead of the control. 
We find that close pairs (r p < 30 fc^j kpc) preferentially contain the 
strongest SFR enhancements (up to a factor of 10 over the control), 
and that galaxies at wider separations show a statistical excess of 
more moderate enhancements (up to a factor of 7 over the con- 
trol). These tests indicate that while the median values are kept at a 
much lower Alog(SFR) due to the relative scarcity of extreme off- 
sets, galaxies in pairs preferentially produce strong SFR enhance- 



ments. However, since this is a test unique to this work, the most 
straightforward comparison to previous studies is a comparison of 
median values. 

The median SFR enhancement in our results is 40% over the 
control; this roughly corresponds to the value of the wide separation 
plateau. The enhancements at smaller r p are between 60% and 80% 
enhancement over the control values. The magnitude of the SFR 
enhancement we find here is broadly consistent with that found in 
previous studies. Ro baina et al.| (2009) find enhancements of or- 
der 80% over the control values for galaxies with r p < 40 fcj] kpc, 
within mass ratios of 1:4, and between redshifts of 0.4 < z < 0.8, 
consis tent with the enhancements seen in our inner peak. |Lin et al.| 
(20071 find that galaxies up to z ~ 1 and r p < 50 kpc are 
enhanced by a factor of 1.9 ± 0.4, relative to a control sample. 
This level of offset is also consistent with the SFR enhancements 
in our sample. |Wong et al.| ^2011) find that the average enhance- 
ment is generally 15-20% for galaxies within 50 fCjl kpc, and 
25-30% for galaxies within 30 hzl kpc for a sample of galaxies 
0.25 < z < 0.75. These values are slightly lower than even our wide 
separation plateau enhancements. However, Wong et al. (20TT) use 
a very generous Av cut of Av< 3000 km s -1 , which will introduce 
a non-negligible fraction of interloping galaxies into their sample, 
their results are a lower limit. Some differences in the observed SFR 
enhancement between works may also be introduced based on dif- 
ferences in the method used to obtain the SFR values, as different 
methods are well known to give different results (e.g., IKennicutt 



& Evans 2012, and references therein). For instance, Wong et al. 
( |2011) use dust-corrected ultraviolet colours as a tracer of SFR, 
while |Lin et al.| {2007) use infrared luminosities to calculate the 
SFRs. Both the effectiveness of the dust corrections and the robust- 
ness of the control sample will influence how robust the enhance- 
ments are to comparisons between studies. With these caveats in 
mind, the relative strength of our results is consistent with those of 
Wong et al. (201 1). It is intriguing that the average SFR enhance- 
ment in galaxy pairs is consistent with our work even at redshifts 
up to z ~ 1 • 

6.4 Metallicity offsets 

There has been very little work on the dependence of metallicity 
on projected separation, and certainly none with the statistics we 
present here. |Kewley et aL| ( |2006| > have a sample of 86 galaxies in 
pairs with r p < 50 JE1 kpc and found that metallicities were only 
significantly diluted at separations closer than 20 fcjjj kpc. Galaxy 
pairs with r,, < 20 A^jj kpc have an average suppression of -0.2 
dex relative to non-pairs. However, Kewley et ah] J2006) uses the 
luminosity-metallicity relation in lieu of a mass-metallicity relation 
to calculate their changes in metallicity, and Elli son et al.| f2008b) 
showed that ~ 50% of this offset was likely due to an increase in 
luminosity due to enhanced SFRs. While our median values are 
not nearly as low as those found in Kewley et al. (2006), (be- 
tween -0.03 and -0.04 dex within 20 hj^ kpc for our total and 
disturbed subsamples respectively), the unbinned points in Fig- 
ure [2] show a significant population of galaxies at offsets of -0.1 
to -0.2 dex, with extreme outliers down to -0.7 dex. Changes 
in metallicity may also be sensitive to internal galaxy parameters 
such as the gas fraction, similar to the dependence for SFR sug- 
gested by |Di Matteo et aL] (2007). Metallicity dilutions have been 
shown to be strongest in low density environments, with metallic- 
ity enhancements present in cluster environments, where the gas 
fraction is likely to be low (Ellison et al. 2009). In support of this 
idea, |Skillman et aT] ( |1996) find that the most gas-deficient galax- 
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ies in the Virgo cluster were also the most metal-rich. Including low 
gas fraction galaxies in the sample may therefore shift the average 
Alog(0/H) to values closer to the control sample. 

Rather than probing metallicity as a function of separation, 
most prior work has has focused on the mass-metallicity relation 
in a sam ple of closer pairs, selecting only galaxies within 25-30 
kpc {Ellison et al.| 



2008b 



Michel-Dansac et al. 2008) where 



the effects of the merger are expected to be the most visible. In 
this context, it has been demonstrated that galaxies in pairs are sys- 
tematically lower in metallicity at fixed mass, relative to a non- 
pair control, and that the magnitude of this effect is approximately 
-0.03-0.05 dex ( |Ellison et al.|2008b||Michel-Dansac et al.|2008| >. 
If we limit our total sample to the galaxies with r p < 30 fcl kpc, 
we obtain a median Alog(0/H) of -0.03 dex; for the disturbed sub- 
sample, the median offset increases to -0.04 dex. |Michel-Dansac| 
et al. ( 2008 1 calculates metallicity offsets in a similar way to this 
work, using the mean metallicity of a set of control galaxies of 
similar mass as the zero point. However, their metallicity offsets 
are only visible in the strongly merging sample, and not the tidally 
disturbed sample. Our division by morphological disturbances does 
not distinguish between these two types, so the higher magnitude 
metallicity offset seen in Michel-Da nsac et al.|p008) (-0.05 dex 
vs. -0.04 dex) is perhaps unsurprising. Therefore, while the form 
of our trend with r p is novel, both the sign and magnitude of the off- 
sets are consistent with previous studies of metallicities in galaxy 
pairs. 

6.5 Simulations 

That the trends we observe in the SDSS data are broadly repro- 
duced by the simple suite of simulations presented here is evidence 
that the merger tracks offered by the models may be an accurate 
framework within which to interpret our results. We have made no 
effort to tune the simulations to our data set, using the same galaxy 
model for all interactions. As mentioned in Section [5] we are not 
reproducing a representative sample of mergers, or fully exploring 
the parameter space of the simulations. For instance, the location of 
the large separation peak seen in Figure [T3fr is partially a function 
of the initial angular momentum given to the galaxy pairs. Fur- 
ther, the simulations in our simple suite do not include any galaxies 
which progress further than ~65 kpc away from their compan- 
ion after the first passage; as our plateau continues until at least 
80 hjg kpc, we expect there to be some fraction of galaxies in our 
sample which are involved in more weakly bound (or unbound) in- 
teractions, which can carry the galaxies to much wider separations 
(Patton et al, in prep). With these caveats in mind, the fact that we 
do reproduce a similar trend to the data, but using a minimum of 
additional assumptions, is very reassuring. This indicates that the 
models do not need to precisely reproduce the exact mergers in the 
data to provide us with a theoretical model of what the galaxies, on 
a statistical level, are doing as they progress through a merger. Us- 
ing these models as an interpretive tool, the increase in the magni- 
tude of offsets in the smallest r p bins is due to the sharp spike seen 
due to galaxies near the end of their merger at those separations, 
when large amounts of metal poor gas are dumped into the nuclear 
region of a galaxy, and a large SFR burst is triggered. This dramatic 
peak due to coalescence at small separations has been seen before 
in both theoretical (e.g., |Di Matteo et al.|2007||Montuori et al.|20T0] 
Torrey et al.|2012) and observational works (e.g., |Larson & Tins- 
ley| 1978| IDonzelli & Pastoriza|1997||Barton et al.|2000||Lambas 
et al. |2003 1| Alonso et al.|2004||2006||Woods & Geller|2007l|Ellison1 
et al.|2008b||20l0l|Darg et al.|2010||Xu et al.|2010[ l. By contrast. 



the wide separation offsets are due to a combination of galaxies 
making their way out after first passage, or possibly back in for a 
second or final passage, depending on the orbits of the interacting 
galaxies. The influence of projected pairs and galaxies that have not 
yet interacted would result in an overall lowering of the magnitude 



7 CONCLUSIONS 

We have used a sample of 1 899 close pairs from the Sloan Digital 
Sky Survey's Data Release 7 to study SFR and metallicity offsets 
as a function of projected separation. We use a simple suite of sim- 
ulations to interpret our results. The main conclusions of this work 
are as follows: 

(i) Galaxies undergoing an interaction show significant metallic- 
ity depressions and SFR enhancements, relative to a control sample 
that is tightly matched in stellar mass, redshift, and local density. 
Galaxy pairs show enhanced SFRs out to projected separations of 
80 hzl kpc, the widest separations in our sample, while significant 
metallicity dilution is observed out to ~60 kpc. Within 30 fc^j 
kpc, SFRs are enhanced, on average, by ~0.21 dex, or by 60%, 
over the control. At separations wider than 30 hz\ kpc, the SFRs 
are enhanced by ~0. 1 dex (25% above the control). Metallicities 
are found to be diluted by -0.02 dex (~5% lower than the control) 
within 60 hzl kpc. 

(ii) When only the morphologically disturbed subsample of 
galaxies is taken, the metallicity offsets are significantly offset from 
the control out to 80 fcj] kpc, and the median offset within 30 hzl 
kpc drops to -0.04 dex (~9% lower than the control). Galaxies 
are visually classified as either visibly disturbed or not visibly dis- 
turbed to attempt to reduce the fraction of galaxies which have not 
yet interacted. 

(iii) The form of the Alog(SFR) and Alog(0/H) offsets as a func- 
tion of r p are primarily driven by galaxies in major mergers. Both 
the more massive and less massive companions in a minor merger 
show similar SFR enhancements to each other and to the major 
mergers. The less massive companions show a minor enhancement 
at all r p , but do not show the same increased SFR enhancement at 
small r p as is visible in the major mergers or more massive compan- 
ions. Furthermore, we find that the most extreme offsets are prefer- 
entially found in major mergers, but that minor mergers are equally 
effective at inducing moderate starbursts in both the less and more 
massive companions. 

(iv) We find that galaxies with strong starbursts, e.g., 10 times 
as strong as the control, are relatively rare, occurring in only 3% 
of galaxy pairs, but are 3 times as likely to be found in the close 
pairs sample (separations < 30 kpc) than in the wide pairs or 
control sample. Extremely diluted metallicities (e.g., 1.8 times as 
metal poor as the control) are rare (3%), but are twice as likely 
to be found in the close pairs sample than in the control. There is 
evidence for synchronous starbursts as a result of a galaxy pair's 
interaction. Galaxies in pairs are at least 17% more likely to show 
significantly enhanced SFR (relative to a scrambled pairs sample) 
if their companion galaxy is also enhanced. 

(v) We use a simple suite of major merger simulations {Torrey| 
|et al.| |2012) to construct a new interpretation of how interactions 
affect the SFRs and metallicities of a galaxy, as a function of pro- 
jected separation. We calculate the Alog(SFR) and Alog(0/H) off- 
sets from the simulations for a qualitative comparison with the 
observational trends. The models are able to reproduce the small 
separation peak in offset value and wider separation offset plateau 
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seen in the SDSS data. The simulations offer the interpretation that 
the wide separation plateau is caused by projection effects blur- 
ring post-pericentre starbursts at wide separations to smaller sep- 
arations. The large offset peak at small separations, by contrast, 
is primarily due to the extreme starburst induced at the end of a 
merger, near coalescence, and is less subject to projection effects. 

Galaxy mergers clearly induce large scale gas inflows in 
galaxies. However, our observations, combined with the theoretical 
framework of hydrodynamical simulations, indicates that these in- 
flows occur on longer timescales than has been previously assumed. 
With metal poor gas reaching the central regions of a galaxy at a 
slower rate, the nuclear gas-phase metallicities will take longer to 
dilute significantly, and the gas reservoir necessary to fuel a signif- 
icant starburst will take longer to accumulate. By the time these 
signatures of gas flow in a galaxy are measurable after the first 
encounter, the two galaxies have progressed out to wide separa- 
tions. As the galaxies reach the end of their merger, the tidal forces 
inducing these gas inflows become significantly stronger, and the 
gas phase metallicities and star formation rates change both more 
rapidly and more dramatically. Future work will quantify the ex- 
tent of separations over which it is possible to see the effects of an 
interaction, and attempt to determine what parameters govern the 
highest SFR galaxies. 
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